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We investigated theoretically electronic and magnetic properties of the perovskite material 
SrCoOs-^ with <5 < 0.15 using a projector-augmented plane-wave method and a Green’s function 
method. This material is known from various experiments to be ferromagnetic with a Curie tem¬ 
perature of 260 K to 305 K and a magnetic moment of 1.5 pe to 3.0 /tb . Applying the magnetic force 
theorem as it is formulated within the Green’s function method, we calculated for SrCoOs-s the 
magnetic exchange interactions and estimated the Curie temperature. Including correlation effects 
by an effective U parameter within the GGA-|-[/ approach and verifying this by hybrid functional 
calculations, we obtained the Curie temperatures in dependence of the oxygen-deficiency close to 
the experimental values. 
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I. INTRODUCTION 

Perovskite materials with the simple structure formula 
ABO 3 attracted attention in the last decades because 
of potential applications in spintronics. Especially, fer¬ 
romagnetic metallic perovskites have rekindled interest, 
since they can be used as electrodes in complex oxide het¬ 
erostructures, while conventional ferromagnetic materials 
such as 3d transition metals are hardly compatible with 
most of the oxides. 

Over the last years, advanced experimental growing 
techniques allow for the growth of heterostructures and 
multilayer systems with a huge variety of properties, e.g. 
multiferroic, magnetoelectric or magnetooptic. Out of 
this class of materials, we concentrated on SrCoOa-a 
(SCO) with a possible oxygen-deficiency 6 . In its pure 
single-crystal composition, the experiments verify ferro¬ 
magnetic behavior up to room temperature and metal¬ 
lic conductance.' Metallic oxides are highly desired for 
contacts and electrodes of the above discussed multi¬ 
functional heterostructures because of a good lattice 
match between the electrode and the top layer oxide. 
Undesired distortions of thin layers could be reduced. Si¬ 
multaneously, it is possible to transport an electric cur¬ 
rent, e.g. in the compounds SrCoi_a;Fea ;03 for x < 0.5, 
which are particularly suitable as electrode materials 
for solid electrolyte oxygen sensors." In addition, com¬ 
plex oxides based on SrCoOa like e.g. Lao. 6 Sro, 4 Co 03 
or SrCoo, 9 Sbo,i 03_5 appear to be particularly useful as 
cathodes for intermediate-to-low temperature solid oxide 
fuel cells. '’d 

In the literature there are several experimental inves¬ 
tigations and theoretical studies on SCO. They all verify 
the ferromagnetic and metallic behavior of this material. 
In theory, in particular, the electronic structure, ’"^ the 
spin state,®’® and possible lattice distortions'® are thor¬ 
oughly discussed. Potze et al.^ show that SrCo 03 ex¬ 
hibits an intermediate spin (IS) state due to the com¬ 


petition of intra-atomic exchange and the cubic crystal 
field. Such an IS state may be understood as a high 
spin state of a Co®d ion antiferromagnetically coupled 
to a ligand hole of Cg symmetry (notation from Ref. 8: 

A more recent study identified this high spin 
state as d®, but it is mixed with several other possible 
spin states.® Also the magnetic moment (ranging from 
2.6 to 3.19is reproduced by theory and close 
to 3 /tb, which agrees very well with the intermediate spin 
state picture {S = 3/2). 

However, there is only one paper, which discussed 
the Curie temperature (Tc) of SrCo 03 from a theo¬ 
retical point of view (Tq ~ 1800K).® This value is 
far too high in comparison to the experimentally mea¬ 
sured values, which vary from 212 K to 305 
This variation might result from different synthesis tech¬ 
niques and defects. Usually, the cubic and stoichiometric 
phase SrCo 03 is prepared from the brownmillerite phase 
SrCo02.5.^ The oxidation is either done by a high oxy¬ 
gen pressure during the heating^ or by electrochemical 
oxidation.'®’'® Furthermore, most samples are polycrys¬ 
talline and the oxygen amount is smaller than the nom¬ 
inal value [5 > 0.05). The cubic structure is only sta¬ 
ble for a narrow range of 5, since SrCoOs-a forms a ho¬ 
mologous series (SrCoO( 3 „_i)/„) for 5 from 0.5 to 0 and 
changes its structure from orthorhombic via cubic and 
tetragonal to cubic again."""'® Even for this cubic range 
close to (5 = 0, the lattice constant increases slightly with 
decreasing oxygen content, *”’'®’'®’'^®’^' whereas the criti¬ 
cal temperature decreases from 292 K to 182 K.'®’^' Any¬ 
way, these values for polycrystalline samples are smaller 
than for a single crystal (Tc = 305 K and 5 = 0.05).’ 
As for the Tc values, the measurements of the magnetic 
moments vary for different experiments between 1.2/tb 
to 2.6AiB-^’“’^^ 

Therefore, we studied in this work the electronic and 
magnetic properties of SrCo 03_5 from ah initio. There 
were strong indications for the IS state proposed by 
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Potze et al.,^ when using in the calculations more ad¬ 
vanced exchange-correlation functionals other than the 
generalized-gradient approximation. Thereby, the theo¬ 
retically estimated critical temperature as well as its re¬ 
duction with increasing amount of oxygen vacancies (Vq) 
agreed well with the experimental observations. 

In the following, we start with a detailed description of 
the studied structure and the applied theoretical meth¬ 
ods. The results are discussed at first for the defect-free 
SCO. Afterwards, we extend our study also to oxygen 
vacancies in SCO. We close this work with a summary. 


II. COMPUTATIONAL DETAILS 

In this work, the calculations were performed for 
a primitive cubic cell of SrCoOa (see Fig. 1). To 
study structural, electronic, and magnetic properties of 
SrCoOa, we performed extensive first-principles calcula¬ 
tions within density functional theory and combined two 
methods, which have been proven to be very reliable in 
this particular context. The projector-augmented plane- 
wave method of the Vienna ab initio simulation pack¬ 
age (VASP)^^“^‘* was used for the calculation of total en¬ 
ergies and structural relaxation, to search for possible 
deviations from the cubic structure. The main results 
were obtained with the Korringa-Kohn-Rostoker Green’s 
function (GF) method HUTSEPOT.“® In both meth¬ 
ods, if not particularly stated otherwise, the generalized- 
gradient approximation (GGA)^*’ was used and we con¬ 
sidered the electronic correlation effects for the d orbitals 
of Co with an additional Hubbard U within the GGA-kG 
approach^' in the implementation of Dudarev et 
The C/eff = U — J was applied as a parameter varying 
from 0 eV to 9 eV. 


A. Structural relaxation 

The plane-wave basis for the VASP calculations was 
taken with a cutoff energy of 460 eV and a P-centered 
(8x8x8) fe-mesh was used. It reproduced the exper¬ 
imental unit cell volume of Ref. 13 (cubic phase; a^ef = 
3.835 A) with a deviation of less than 3 % for the studied 
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FIG. 2. Experimental lattice volume as a function of oxygen 
concentration from several references. The change of the vol¬ 
ume is AV = {V — a?ef)/fl?ef with respect to Oref = 3.835 A 
(black dashed line). The nominal and net values of Ref. 15 
represent the provided and measured oxygen content. 


range of UeS- In the same order of magnitude or smaller 
are the differences (with respect to Oref) to the single¬ 
crystal lattice constant' Osngi = 3.8289 A (AV = —0.47% 
with V = (osngi)^) or possible variations due to oxygen 
vacancies observed e.g. in Ref. 15 (AV = -1-2.2%) or 
Ref. 16 (AV = -1-0.56 %) (see overview in Fig. 2). For the 
sake of completeness, we compared our GGA results with 
hybrid functional calculations using the common HSE03 
functional,^® available within VASP but at the moment 
not available in the GF method. In the HSE03 func¬ 
tional, the mixing a = 0.25 and the screening parameter 
/i = 0.3 A“^ are chosen for the Hartree-Fock exchange 
energy. Because of the high computational demands of 
such type of calculations the cutoff energy was reduced to 
400 eV and the F-centered fc-mesh was only 4x4x4. The 
lattice relaxation provided for a cubic phase the equilib¬ 
rium lattice parameter ohseos = 3.833 A, which is in a 
very good agreement with the experimental values shown 
above. Although we did not observe any indication of 
a tetragonal distortion in our GGA-I-G calculations, we 
found by applying the HSE03 exchange functional that 
the a/c ratio deviated by few percent from one. Since the 
experiments for stoichiometric SrGoOg did not observed 
a tetragonal unit cell and the electronic structure did not 
change substantially with respect to the cubic phase, we 
stayed in the main part of the paper consistently with the 
cubic SCO and the lattice constant Oref = 3.835 A and 
investigated a hydrostatic volume variation separately in 
Sec. BID. 


B. Magnetic properties with the Green’s function 
and Monte Carlo method 




3 


The calculations with the GF method were performed 
within the full charge density approximation, which takes 
into account the non-sphericity of the charge density and 
improved the accuracy of calculations for complex unit 
cell geometries. 

The GF method, which is based on the multiple scat¬ 
tering theory, allows the calculation of the magnetic ex¬ 
change interactions Jij between the magnetic atoms at 
site i or j by using the magnetic force theorem.'^'' Therein, 
the magnetic moments at i and j of an ordered magnetic 
structure are rotated against each other about a small 
angle. From the resulting energy variation follows 

1 r 

Jij = deImTri(A,r*^AjT;J* -h A,t*^AjT^*) , (1) 

where is the scattering path operator and A^ = 
— i~^ the difference between the spin-dependent sin¬ 
gle scattering operators of site i. The trace is taken 
over all relevant angular momentum quantum numbers 
L = {I, m). 

The Jij from (1) enter a classical Heisenberg Hamilto¬ 
nian 

H = — Jij Si ■ Sj , ( 2 ) 

i,j 

with the magnetic moments Si and Sj. The model 
Hamiltonian (2) allows an estimate of the critical tem¬ 
perature T with a Monte Carlo (MG) simulation. Here, 
2 ^mc mark the transition to any kind of magnetic 

ground state, which might result from the calculations. 
We cross-checked our results with the mean-field approx¬ 
imation (MFA) and the well established random phase 
approximation (RPA).^'^ All three methods showed a sim¬ 
ilar tendency in the variation of T with respect to the 
considered influences like electron correlations or oxygen 
vacancies. Since the MFA is known to overestimate T, we 
restricted the discussion to the MG results. For the par¬ 
ticular MC simulations, we choose 8000 magnetic atoms 
in a cluster with periodic boundary conditions. The tem¬ 
perature was reduced in steps of 5 K starting from a high- 
temperature disordered state above . At every tem¬ 
perature T, we assumed that the thermal equilibrium was 
reached after 20 000 MC steps and after 20 000 additional 
MC steps thermal averages were calculated. was 

then obtained from the htting of the temperature depen¬ 
dency of the magnetic susceptibility, cross-checked by the 
temperature dependence of the saturation magnetization 
and the heat capacity. The obtained transition temper¬ 
atures were determined within a numerical uncertainty 
range of ±5 K. Further computational details of our MC 
scheme can be found in Refs. [32-34]. From the orienta¬ 
tion of the magnetic moments at low temperatures and 
the spin-spin correlation function, we deduced the mag¬ 
netic ground state of the simulated system, which was 
not only ferromagnetic (FM) with the Curie temperature, 
Tq, as observed from the experiments, but also antifer¬ 
romagnetic (AFM) with the Neel temperature, Tn, or 


a more complicated non-collinear ferri-magnetic spin ar¬ 
rangement (FiM) with a transition temperature, T, and 
potentially with a saturation magnetization of zero. We 
introduced those labels for the different magnetic tran¬ 
sitions for the sake of clarity throughout our work and 
speak of T if the ground state is not clarified. 


In our calculations for a ferromagnetic ground state 
(reference state - RSfm), we found besides the magnetic 
cobalt ions small induced moments at the oxygen sites. 
Since those moments disappear usually at and above 
the magnetic transition temperature, their magnetic cou¬ 
pling can lead to a wrong estimation of the critical tem¬ 
perature. We compared the results of (1) for the FM ref¬ 
erence state RSpM with the magnetic coupling constants 
obtained at a high temperature paramagnetic (PM) refer¬ 
ence state (RSpm)- Such state can be modeled success¬ 
fully with the disordered local moment (DLM) theory 
using the GF method.'*®’^® Here, arrangements of local 
magnetic moments {Si} at the sites i are thought to fluc¬ 
tuate independently. Above T, the orientations of those 
local moments are randomly distributed and the average 
magnetization per site is zero. Hence, the induced mo¬ 
ments vanish. From the computational point of view, the 
coherent potential approximation (CPA)'^^ as it is imple¬ 
mented within the multiple scattering theory** can be 
used to model the susceptibility or the electronic struc¬ 
ture for such disordered magnetic systems.*®’**’ With re¬ 
spect to the calculation of the Jij in (1), the scattering 
path operators f*-’ of the perfect magnetic ordering will 
be exchanged with the ones of the effective GPA medium 
Tq . The random orientation of the magnetic moments 
in the DLM picture causes usually an increase in the 
size of the orbitals. So, the magnetic coupling is usually 
stronger in case of a DLM calculation since the strength 
of the magnetic coupling is proportional to the overlap 
between the contributing orbitals. This means for the 
resulting transition temperature that it can become ei¬ 
ther larger or smaller depending on the type of magnetic 
exchange, e.g. a reduction for antiferromagnetic (super-) 
exchange. 


To describe the oxygen-deficient SrCoOa-^ in the GF 
method, we applied also the GPA using a certain amount 
(i5) of empty spheres at the oxygen sites to mimic Vq. 
However, vacancies may lead to substantial relaxations 
of the underlying crystal structure. We did not account 
for such structural deformations in our CPA calculations 
but investigated their impact on the magnetic interaction 
using a supercell approach with the GF method. We 
found only minor changes in the exchange constant values 
and, therefore, the discussion below reports results from 
the ideal cubic structure. 
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FIG. 3. (Color online) Atomic and spin-resolved density of 
states of SrCoOa with contributions of each species and spin- 
up (upper panel) and spin-down (lower panel) obtained with 
the GF method (GGA). The oxygen p states at all three oxy¬ 
gen ions in the unit cell are three fold degenerate (only one is 
shown). Main feature of the LDOS are similar in VASP (not 
shown). 


III. RESULTS 

A. Defect-free SCO 

We investigated at first the electronic properties of 
SrCoOa for the experimental lattice constant. Due to 
the oxygen octahedron, which surrounds the Co ion, the 
d states of cobalt experience a crystal field splitting of 
cubic symmetry, which results in three t 2 g and two Cg 
degenerated states. The orbitals corresponding to the Cg 
states are oriented along the coordinate axes pointing to 
the oxygen ions and those corresponding to the t 2 g states 
are pointing to the next nearest neighboring Co ions (see 
Fig. 1). 

The local density of states (LDOS) within the GGA 
shows for VASP and HUTSEPOT an almost fully oc¬ 
cupied majority spin channel and a pronounced peak of 
the Co t 2 g spin-down states at the Fermi energy Ep (see 
Fig. 3). All other cobalt 3d states, ej, and e|, are 
below the Ep and smeared over a large energy range due 
to a strong hybridization with the p states of oxygen. 
This is in a good agreement with previous results*^ but 
contradicts the IS state picture with particular occupied 
and unoccupied t 2 g Co states (t^gCg)-^ 

The total moment of this IS state model would be the¬ 
oretically 3 ^b (s = 3/2). Although the total magnetic 
moment of /i = 2.281 /tb calculated with the GF method 
and the GGA functional was smaller than in the IS state 
model, it agreed well with the range of the experimen¬ 
tally observed magnetic moments (p, « 2 /tb to 2.5 pp). 
The main contribution to the total moment originated 
correctly from the Go ions with only small induced mo¬ 
ments of 0.15/iB and 0.04 ^b at the oxygen and stron¬ 
tium ions, respectively. The discrepancy between the IS 
state model and the experiments was always attributed 



FIG. 4. (Color online) Spin-resolved density of states of 
SrCoOs in the spin-up (upper panel) and spin-down chan¬ 
nels (lower panel) obtained with VASP (HSE03 functional) 
and GF method (GGA-tC with U = 5eV). Only the non 
degenerated Co 3d states are shown. 


to possible defects. However, when the IS state model is 
valid, theoretical calculations of a defect-free SGO should 
reproduce the total magnetic moment. 

Possible shortcomings might originate from the GGA 
electron correlation functional, which often lacks a suf¬ 
ficient description of localized transition metal d states. 
For comparison, we calculated the total magnetic mo¬ 
ment and the DOS of SCO also with the hybrid functional 
HSE03 available in VASP. The obtained total magnetic 
moment of 2.9 pp within VASP and the HSE03 functional 
agreed well with the IS state model. Additionally, the 
DOS showed also the expected orbital occupation t^gCg of 
the Co states: two unoccupied t^g {dxz and dyz above Ep) 
and one occupied state {dxy below Ep) (see Fig. 4). 

We have to note that due to the symmetric cubic struc¬ 
ture, the particular localization of is arbitrary and 
depends on the starting point of the self-consistent cal¬ 
culation. Another localized Co d^ state is also possible 
and was observed during the calculations. However, for a 
consistent description we continued throughout this work 
with one particular configuration (the singlet d)^y state). 

In order to calculate the critical temperature, the 
magnetic exchange interactions for the nearest neighbor 
atoms were calculated with the GF method and the GGA 
functional. For RSfm the magnetic interaction param¬ 
eters show mainly a ferromagnetic (positive) coupling: 
strong between two adjacent Co ions and much weaker 
between the Co ions and the induced magnetic moments 
of the surrounding oxygen ions (see Fig. 5). The coupling 
to the Sr ions was one order of magnitude smaller and was 
ignored in the following discussions. After the 8th shell 
(d = 7.7 A), most of the coupling constants decay fast 
while other long-range interactions reflect the metallic 
character of SrCoOa. Along the straight exchange paths 
(along a particular coordination axis, also marked with 
an asterisk in Fig. 5) the coupling via the oxygen ions 
remains stronger. The same tendency can be observed 
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FIG. 5. (Color online) (a) Schematic view of the orientation 
of the first magnetic exchange interactions Jij only between 
the Co ions. The non magnetic Sr ions are not shown. The 
notation Joi, J 02 , ... means nearest neighbors (NN), next 
NN, etc. Below, the values for the coupling between Co-Co 
and Co-O are given in (b) and (c), respectively. Both calcu¬ 
lated with the GF method (GCA) for a FM and PM (DLM 
model) reference state, respectively. The asterisk marks those 
coupling constants which strictly follow only one direction in 
space (along a single cube edge in (a)). The abscissa in (b) is 
non-continuous, due to the large differences in the values. 


for the coupling between the Co ions in the DLM model 
{Jij calculated with RSpm)- Only the nearest neighbor 
interactions became larger while the other coupling con¬ 
stants were reduced. At the oxygen ions, the induced 
magnetic moments vanished in the DLM model and no 
exchange interactions were found between them. In the 
Monte Carlo method, we took into account all coupling 
constants up to the distance of 15.34 A. The resulting 
Curie temperature did not agree with the experimental 
results (within the DLM model = 771K compared 
to « 280 K). The reason for this high might be 
an overestimation of the magnetic coupling between the 
Co atoms. The increase in the localization of 3d states by 
using electron correlation corrections or more advanced 
functionals should decrease as well the orbital overlap 
and thereby reduce the exchange coupling. 

In summary, the DOS, the total magnetic moment and 
the critical temperature, illustrate the need of electron 
correlation corrections for the description of the proper- 
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FIG. 6. (Color online) Contour plots of the local density 
of states (LDOS) in the ground state for the Co d states in 
dependence of the correlation parameter 17eff (ordinate) ob¬ 
tained with the GF method. They are spin-resolved (left: up, 
right: down) and collected according to their initial degener¬ 
acy. For higher UeB (not shown), the energy shift continues 
and there are no significant changes to observe. 


ties of SCO. Unfortunately, the calculations with hybrid 
functionals consume large amounts of computational re¬ 
sources and were also at the moment not available within 
HUTSEPOT, which was needed for the calculations of 
the Jij. However, the parametric GGA-I-C/ approach 
was less time consuming and provided similar results, 
as it is discussed in the following Sec. IIIB. For that 
reason, we stayed in the remaining of this work with 
the GGA-hU method. Additionally, we describe also 
the oxygen-deficient material, where we expect a reduc¬ 
tion of the critical temperature and a slightly enhanced 
volume.^'’ 


B. Electronic correlation corrections 

The common way to consider the electronic correla¬ 
tions within the GGA-|-f7 method is to optimize the value 
of the repulsive Ugs with respect to the experimental data 
for structural and magnetic properties of the system. 
As we see from former studies, the value of UeS might 
range from 2.5 eV (Ref. 10) to 8eV (Ref. 7). The con¬ 
strained random-phase approximation provided a value 
of t/ = 10.83 eV and J = 0.76 eV for the Co d states'’ 
which seems to be too high for a metallic system. 

Since the correct value of UeS is hardly to estimate 
from first-principles and the above reference values scat¬ 
ter quite a bit, we investigated the electronic structure 
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and the occupation of the Co d states for the whole range 
of UeB from 0 eV to 9 eV with the GF method (see Fig. 6 ). 
Interestingly, for the first few steps of the calculations 
{UeS < 1.5 eV) the d states preserve their degeneracy in 
t 2 g and Cg states. At the Fermi energy Ep, the large peak 
of the Co t 2 g states does not move due to an interplay 
of the Coulomb exchange and the crystal-field energy. 
Only for a larger UeS, the degeneracy is lifted as it was 
observed before in the DOS calculated with the HSE03 
(see Fig. 4). A singlet state (d^y) becomes occupied while 
the doublet {d^^ and is pushed above Ep. On the 
other hand, in the spin-up channel the orbitals remain de¬ 
generate for the whole range of Ues and become strongly 
localized (see much higher contrast for d^y in Fig. 6 ). It 
matches well with the IS state model. 

A similar loss of degeneracy in the 3d states was also 
observed in calculations with VASP and GGA-bC, but at 
higher values for Ues than in the GF method (happening 
between 5eV to 7eV, not shown). Furthermore, we ob¬ 
tained a very good agreement of the electronic structure 
around the Ep, calculated either with the GF method 
at UeS = 5eV or with VASP (HSE03) (see Fig. 4). For 
states further below Ep, the differences become larger 
but those states contribute only little to the orbital over¬ 
lap needed for the calculation of magnetic exchange pa¬ 
rameters with the magnetic force theorem. 

For the understanding of the critical temperature, we 
calculated with the GF method the magnetic coupling 
parameters Jij for the FM and PM (DEM model) ref¬ 
erence state, respectively (see Fig. 7). Both reflect the 
change in the degeneracy of the Go d state in the LDOS. 
Due to the splitting of the t^g states into the degenerated 
doublet (d);^ and d^^) and a singlet (d;^:^^), the Jy with ei¬ 
ther (II) in a; and y direction or (_L) in z direction were dif¬ 
ferent for UeS > 1.5 eV. For the ground state calculation 
(RSfm), the most dominant coupling constants, 
and are strongly ferromagnetic (see left hand side 

in Fig. 7). Both interactions are mediated by oxygen ions 
between Co ions forming either a Co-O-Co or a Co-O-Co- 
0-Co chain (see inset in Fig. 5). Those bonds connect 
mostly the O p states with the Co Cg states (cr bonds). 
This typically antiferromagnetic (AFM) superexchange 
is suppressed by the metallic character of SrCoOa and 
we observe band magnetism. In the degenerated param¬ 
eter region, the coupling becomes even stronger with Ues 
until 1.5 eV because the Co e| and states are either 
pushed below or above Ep, respectively. This increases 
the exchange splitting and, therefore, the magnetic cou¬ 
pling. In contrast, the coupling between two Co ions 
enclosing a 90° angle (d^°“^°) is small and AFM, while 
the next coupling (d^°“^°) is very weak compared to the 
other interactions. 

After that, in the symmetry broken regime, the com¬ 
peting superexchange overcomes the band magnetism 
and we observed an increasing localization of the d states 
(see stronger contrast in Fig. 6 ). This reduces, in general, 
the overlap of the orbitals and the magnetic exchange in¬ 
teractions. It is visible e.g. for Jqi and J 04 in Fig. 7 but 


RSfm 


RSpM 



FIG. 7. (Color online) The magnetic exchange interactions 
jCo-Co function of UeS for a ferromagnetic (left) or para¬ 
magnetic (right) reference system. They are subdivided into 
two groups, II aud T, for only in the x-y-plane and 

those with contributions also in x direction (directions given 
in Fig. 5(a)). The 7^°'^° (not shown) does not split due to 
f/eff and are always small (< 0.3meV). 


also for other Jqj (not shown). On the other hand, the 
modifications in the coupling constants are much more 
complex due to the changing LDOS. However, some sim¬ 
ple tendencies can be observed, e.g. due to the local¬ 
ization of the Go dxy states, the magnetic coupling for 
UeS < 4eV in the x-y-plane (|| contributions in Jqi) be¬ 
comes smaller than the out-of-plane (T) contributions. 
The T parts of even change their character from 

AFM to FM. Another significant change for all coupling 
constants is visible at UeS = 4.5 eV and might be corre¬ 
lated with vanished states in the dxy LDOS indicated by 
a loss of contrast (see Fig. 6 ). At the end of the shown 
range, the strength of the nearest neighbor magnetic ex¬ 
change coupling was only half of its starting value. It re¬ 
duces further, for even higher UeS (not shown) and leads 
to an undesired antiferromagnetic ground state. 

The overall tendencies for the RSfm calculations were 
in general also observed in the DLM picture (RSpm), al¬ 
though the changes were much stronger, e.g. at Ues = 
6 eV the Jqi were reduced to zero and J 04 is strongly 
AFM. Furthermore, the loss of degeneracy was visible 
already for smaller UeS (see Fig. 7). Both changes are 
explained by the larger extent of the Go d orbital due to 
the random distribution of the magnetic moments in the 




















7 



FIG. 8. (Color online) The critical temperatures for the mag¬ 
netic transition (left axis) and the total magnetic moment per 
functional unit at the ground state (right axis) as a function 
of ?7eff. The Monte Carlo calculations show a change in the 
magnetic ground state from clear ferromagnetic behavior to 
a FiM situation at UeB ~ 2.1 eV and 4eV for the RSpm and 
RSfm calculations, respectively (see text). The gray shaded 
area indicates the experimental Tc range. 


DLM theory as already stated in section II. The increas¬ 
ing orbital overlap enhances, on the one hand the AFM 
superexchange and on the other hand, alters the com¬ 
petition between the crystal field and Coulomb energy 
which restored the degeneracy for small UeS- 

The critical temperatures obtained from the Monte 
Carlo simulation using the magnetic exchange parame¬ 
ters in the RSfm or with the DLM model follow in gen¬ 
eral a similar tendency as the nearest neighbor coupling 
constant (see Fig. 8). They show for both sets 

of magnetic coupling parameters a linear increase up to 
[/eff = 1.5 eV. The critical temperatures obtained with 
the JojXRSfm) remained around 750 K with increasing 
UeS up to « 3eV, and drop down sharply in the follow¬ 
ing while the DLM results decrease linearly immediately 
above 1.5 eV. So, T^^{Ues) calculated for RSpm reaches 
the experimentally relevant range already for a smaller 
?7eff = 2.5eV than for RSfm (4.5eV). A reason for the 
smaller UeS was already discussed above for the Jij - the 
larger overlap of the orbitals in the more realistic DLM 
model. Additionally, the RSpm calculation showed af¬ 
ter the kink an almost constant T^^{Ues) = 240 K for 
a larger range of UeS parameters (2.75 eV to 4eV). On 
the other side, the ground state calculation at RSfm with 
the GF method (GGA-I-C) returns for UeS » 5eV a 
inside the experimental range. This observation matches 
well with the above discussion about the electronic struc¬ 
ture and the correspondence of GF method (GGA-|-?7) 
and VASP (HSE03) calculations at Ues ~ 5eV (see 
Fig. 4). 

Furthermore, we account in Fig. 8 for the magnetic 
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FIG. 9. (Color online) LDOS of the oxygen along the 2 di¬ 
rection (Oz) for GCA, UeB = 2.5 eV and with 5at.% oxygen 
vacancies (Vo). 


ground state observed in the Monte Garlo study with 
different symbols, having either a FM or FiM ground 
state. The DLM theory predicted only a FiM ordered 
ground state for the range of ([/es « 2.5 eV) equiv¬ 
alent to experimental results. This is, however, still phys¬ 
ically reasonable. While at the critical temperature the 
induced moments at oxygen might be zero, they will ap¬ 
pear at lower temperatures. This leads to a ferromagnetic 
ground state obtained with the Jij (RSfm) for even larger 
UeS (see blue curve in Fig. 8). 

On the other hand, the total magnetic moment in¬ 
creases monotonously in the whole UeS range with dif¬ 
ferent linear slopes (see red line in Fig. 8). Only for 
changing degeneracy for Ues > 1.5 eV the total moment 
decreases slightly due to a reduction of the induced mo¬ 
ments to 0.11/ip for the two oxygen ions (Ox and Oy), 
which lie in the same cc-y-plane as the Co ions. In the fol¬ 
lowing, the linear slope changed around Ues ~ 2.5 eV and 
UeB ~ 4.5 eV, which correspond to disappearing peaks in 
the LDOS (see Fig. 6). At the latter UeS, the total mo¬ 
ment is 3 fj,B and matches well on the one hand side with 
the IS state model and on the other hand with the HSE03 
calculation in VASP. 

We found in summary that the PM model (DLM the¬ 
ory) yields for the critical temperature a good qualitative 
agreement with the measurements *4'’ only for a small 
correlation parameter of Ues = 2.50 eV to 2.75 eV. Still, 
a exact comparison with the measured Tq remained com¬ 
plicated due to the different experimental setups, single¬ 
crystalline or poly crystalline samples, varying growing 
techniques or different oxygen content. 


C. Effects of oxygen vacancies 

We used the CPA of the GF method to substitute 
the oxygen sites with a certain concentration of empty 
spheres, modeling the oxygen vacancies. For a low oxy¬ 
gen deficiency in SrGoOa-i up to 5at.%, there is no 
experimental evidence for an ordering of Vq. So, we 
assumed randomly distributed oxygen vacancies which 
matches well with the concept of the GPA. 
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FIG. 10. (Color online) The theoretical critical tempera¬ 
ture (DLM model) in dependence of the oxygen content in 
SrCoOs-i compared with experimental references. The GGA 
calculations shown as dashed blue line and GGA-|-t/ calcula¬ 
tions shown as differently colored solid lines. The open and 
full circles in the corresponding colors indicate the calculated 
data points and the magnetic ground state as shown in Fig. 8. 
An additional antiferromagnetic state (AFM) is marked with 
a half filled circle (see text). 


The electronic structure of the oxygen ions is mainly 
dominated by the strong hybridization with the Co d 
states (see Fig. 3). By including of few atomic percent of 
oxygen vacancies (Sat. %) in SCO, the unoccupied peak 
above Ep is shifted to higher energies (see Fig. 9), a sim¬ 
ilar effect as increasing UeS (see GGA results for com¬ 
parison in Fig. 9). Both lead to an enhancing orbital lo¬ 
calization - either removing effectively oxygen from the 
lattice or reducing the electron hopping between the Co 
ions. Anyway, these unoccnpied oxygen states can be 
interpreted as the ligand hole of eg symmetry, which is 
expected for the IS state ^ 2 g 6 g-* I* is stabilized by the 
ligand hole state d^L^^ where the hole couples antiferro- 
magnetically to another e| in d®: t^gCg- We found this 
configuration also in our LDOS calculation with almost 
fully occupied orbitals for d^ and t^g (see 6 ). Although 
our method is not directly comparable with the dynamic 
mean-held theory (DMFT) method of Kunes et alP , we 
note that the same spin conhguration appeared also in 
their calculation with the highest multiplet weight. 

We calculated again the magnetic transition tempera¬ 
ture T^^{S) based on the ab initio magnetic exchange 
coupling constants in RSpm with the MC method but 
included different amounts of oxygen vacancies 5. Al¬ 
though we observed at the relevant critical tempera¬ 
tures no ferromagnetic ordering, calculations with the 
Jij in RSpM revealed a ferromagnetic ground state at 
low temperatures for all considered vacancy amounts (not 


shown). The results were compared with the experimen¬ 
tal measurements.^’^® The experiments with polycrys¬ 
talline samples^® varied the oxygen content by ± 6 at.% 
with respect to the stoichiometric sample and found de¬ 
creasing Tc with reduced oxygen amount (see Fig. 10). 
A positive value means excess oxygen, which might oc¬ 
cupy unknown interstitial sites or form more complicated 
point defects. This is beyond our aims and we restricted 
the current study to model only oxygen deficiency by 
the introduction of oxygen vacancies into the cubic unit 
cell. In contrast, the single-crystalline Tq was found to 
be remarkably higher, but the lack of more data points 
complicates the comparison.^ 

In theory, we expected also a reduction of 
since we found that the localization of the orbitals 
through oxygen vacancies was like applying an UeB pa¬ 
rameter. Therefore, the magnetic interactions in the 
Go-O-Co bonds were weakened as well. This expecta¬ 
tion was generally fulfilled for the calculated T^^{S) 
at different Ues' the GGA exchange functional alone, 
UeB = 2 eV, 2.25 eV, 2.5 eV or 2.75 eV (see Fig. 10). They 
all show an increase of the critical temperatures by in¬ 
creasing the oxygen content from 6 = 0.3 (lOat.%) to¬ 
wards stoichiometric SCO while their slopes varies qual¬ 
itatively between UeS = 2.25 eV and 2.5 eV. The slope 
remained equal until Uqb = 2.25 eV and even partially 
2.5 eV but became for larger UeB reduced and matched 
well the trend of the experimental results. We note that 
for UeB = 2.5 eV only the combined effect of UeB and 
a larger amount of Vq is enough to obtain the experi¬ 
mental trend, whereas UeS = 2.75 eV lies in the region 
of constant T“*^([/eff). For the latter, the influence of 
the oxygen vacancies alone to T^^{6) was visible and it 
agreed well with the measurements (see Fig. 8 ). 

However, these results were just valid in a small range 
of oxygen-deficiency, otherwise SrCoOa-a becomes un¬ 
stable and forms different structure. This is also visi¬ 
ble in the observed magnetic ground states depicted in 
Fig. 8 . At 15at.% oxygen-deficiency and UeS > 2.5 eV, 
the resulting magnetic ground state was antiferromag¬ 
netic. This chemical composition is already close to 
the ordered brownmillerite structure SrCo 02 . 5 , which 
is antiferromagnetic with a high Neel temperature of 
Tn = 570 K.'^ Although the cubic structure is not the 
appropriate equilibrium structure at this oxygen concen¬ 
tration, the difference in chemical composition might lead 
already to an AFM order with a higher critical temper¬ 
ature. Such a transition from a FM to an AFM ground 
state can be caused in SCO e.g. by strain.^® A strain in 
general varies the unit cell volume and might therefore 
modify the magnetic coupling as well. 

D. Hydrostatic volume changes 

Up to now all calculations were performed with the 
fixed lattice constant Oref, but measured lattice constants 
indicate a volume expansion in consequence of oxygen- 
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FIG. 11. (Color online) The theoretical critical temperature 
(DLM model, GGA+t/) of SrCoOa-i in dependence of the 
volume change AV (filled squares; lower axis; reversed to 
match other figures). The volume change is correlated with 
the oxygen-deficiency (upper axis) via the experimental re¬ 
sults of Ref. 15 (see text). Both effects are also combined 
(open squares). Ues is 2.5eV (blue) or 2.75eV (black), re¬ 
spectively. The magnetic ground state is always FiM. 

deficiency (see examples in Fig. 2). To estimate the influ¬ 
ence of this volume enhancement to the critical tempera¬ 
ture f“C(Ay), we choose from Fig. 2 the largest volume 
expansion as reference.''”’ So, we scaled the lattice param¬ 
eters of the cubic unit cell for defect-free SCO up to 1.8 % 
and calculated in RSpM (see filled squares and the 
lower axis in Fig. 11). For UeS = 2.5eV, the increasing 
distance between the Co atoms reduces their magnetic 
interaction and the critical temperature, again similar as 
Ueff or Vo- Hence, T''’^®(AV) remains for UeS = 2.75eV 
almost constant. 

On the other hand, the volume expansion is corre¬ 
lated to a particular oxygen content (see Fig. 2). For 
simplicity, it is linearly interpolated from Ref. 15 to 
derive in Fig. 11 the upper axis (open squares). The 
combination of volume expansion and oxygen-deficiency 
leads for UeS = 2.5 eV to a similar quantitative curve as 


for T^^(5) but the observed kink appears already for 
5 = 0.05 or AV = 0.6%, respectively. In contrast, the 
qualitative and quantitative agreement of the variation 
of for UeB = 2.75 eV is still in place. Finally, we 
note that although the good agreement with Ref. 15 of 
the theoretically obtained , the variation AV was 
the upper boundary of Fig. 2. For example, in Ref. 18 
the oxygen-deficiency of 5 = 0.16 with respect to their 
value at 5 = 0 is correlated to a much smaller volume ex¬ 
pansion of AV = 0.6%. As a result the slope of T^^(5) 
changes drastically as well. 

IV. CONCLUSIONS 

We conclude that for cubic SrCoOa-^ the inclusion of 
temperature effects via the DLM model and a small cor¬ 
relation correction parameter of Ueg « 2.75 eV is neces¬ 
sary to describe the measured magnetic properties, such 
as the magnetic moment or, in particular, the Curie tem¬ 
perature. Such values of UeB are expected for a metallic 
system. In the studied compound, we observe mainly the 
band magnetism, which is reduced by oxygen-mediated 
superexchange. Furthermore, our calculations agree with 
the atomic multiplet calculations and the picture of the 
intermediate spin state of Potze et al.^ 

On the other side, oxygen vacancies can drastically 
alter the magnetic properties.The simple means of 
the coherent potential approximation were successful 
to model qualitatively and quantitatively the experi¬ 
mentally observed reduction of T induced by oxygen 
vacancies.'® They are one of the most important types of 
defects in oxides - even single crystalline samples might 
be not completely stoichiometric.' The randomly dis¬ 
tributed vacancies weaken the Co-O-Co bonds and re¬ 
duce the exchange coupling, similar to electronic correla¬ 
tions. The same behavior was also observed for a possible 
increase of the unit cell volume. 
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